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Abstract 

We present a physically reasonable source for an static, axially- 
symmetric solution to the Einstein equations. Arguments are pro- 
vided, supporting our belief that the exterior space-time produced by 
such source, describing a quadrupole correction to the Schwarzschild 
metric, is particularly suitable (among known solutions of the Weyl 
family) for discussing the properties of quasi-spherical gravitational 
fields. 
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1 Introduction 



As is well known, Weyl solutions [1] represent the family of all static and 
axially-symmetric exterior solutions to the Einstein equations. Since there 
are as many different Weyl solutions as there are different harmonic functions 
(see next section), then the obvious question arises: what is the exact vacuum 
solution to the Einstein equations corresponding to a given non-spherical, 
static axially symmetric source (an ellipsoid, say)?. 

If the field is not particularly intense (r 2M) and the deviation from 
spherical symmetry is sligth, then there is not problem in representing the 
corresponding field (both inside and outside the source) as a suitable pertur- 
bation of the spherically symmetric exact solution. However, as the object 
becomes more and more compact, such perturbative scheme will eventually 
fail. 

Indeed, as it is well known [2], the only static and asymptotically-flat vac- 
uum space-time possessing a regular horizon is the Schwarzschild solution. 
For all the others Weyl exterior solutions [1], the physical components of the 
Riemann tensor exhibit singularities at r = 2M. Therefore, it is intuitively 
clear that as r approaches 2M and gravitational field becomes stronger, the 
properties of sources of Weyl space-time should start drastically to differ 
from the properties of spherical sources [3]. It is important to keep in mind 

2 



that this sharp difference in the behaviour of both types of sources (for very 
high gravitational fields) is independent on how small, multipole moments 
(higher than monopole) of the Weyl source, are. This is so because, as r 
approaches 2M, any finite perturbation of the Schwarzschild space-time be- 
comes fundamentally different from any Weyl solution, even when the latter 
is characterized by parameters whose values are arbitrarily close to those cor- 
responding to the spherical symmetry This point has been stressed before 
[4], but usually it has been overlooked. 

Furthermore in a recent work [5] it was shown that for a non-spherical 
source (even in the case of slight deviations from spherical symmetry), the 
speed of entering the collapse regime decreases substantially, as compared 
with the exactly spherically symmetric case. In the same order of ideas it 
has been shown [6] that small departures from sphericity, produce significant 
decreasing (increasing) in the values of active gravitational mass of collaps- 
ing (expanding) spheres, with respect to its value in equilibrium, enhancing 
thereby the stability of the system. Also, the sensitivity of the trajectories 
of test particles in the 7 spacetime, to small changes of 7, for orbits close to 
2M, has been brought out [7] . It is important to stress that all these effects 
take place for strong gravitational fields, but for r > 2M. 

These works [5], [6], [7] were done using as Weyl solution, the so-called 
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gamma metric (7-metric) [8]. This metric, which is also known as Zipoy- 
Vorhees metric, belongs to the family of Weyl's solutions, and is continuously 
linked to the Schwarzschild space-time through one of its parameters. The 
motivation for this choice was that the exterior 7-metric corresponds to a 
solution of the Laplace equation (in cylindrical coordinates) with the same 
singularity structure as the Schwarzschild solution (a line segment [8]). In this 
sense the 7-metric appears as the "natural" generalization of Schwarzschild 
space-time to the axisymmetric case. 

However, the existence of so many different (physically distinguishable 
[9]) Weyl solutions gives rise to the question: which among Weyl solutions is 
better entitled to describe small deviations from spherical symmetry?. 

Although it should be obvious that such a question does not have a unique 
answer (there is an infinite number of ways of being non-spherical, so to 
speak), we shall invoke a very simple criterion, emerging from Newtonian 
gravity, in order to choose our solution. 

Indeed, in order to answer the question above, let us recall [10] that most 
known Weyl solutions, present a drawback when describing quasi-spherical 
space-times. It consists of the fact that its multipole structure is such that 
multipole moments, higher than quadrupole, are of the same order as the 
quadrupole moment. Instead, as it is intuitively clear, the relevance of such 
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multipole moments should decrease as we move from lower to higher mo- 
ments, the quadrupole moment being the most relevant for a small departure 
from sphericity. Thus for example in Newtonian gravity, multipole moments 
of an ellipsoid of rotation, with homogeneous density, and axes (a,a,b) read: 

D 2n = (-2) n 3Ma 2n e n (l -e/2) n /[(2n+ l)(2n + 3)], e = (a-b)/a, 
D 2n+ i = 

because of the factor e n , this equation clearly exhibits the progressive de- 
creasing of the relevance of multipole moments as n increases. 

Thus, in order to describe small departures from sphericity, by means of 
exact solution to the Einstein equations, we would require an exact solution 
whose multipole structure shares the property mentioned above. Fortunately 
enough, such solution exists [10]. Indeed, there is one (exact) solution of the 
Weyl family, which may be interpreted as a quadrupole correction to the 
Schwarzschild space-time (see below). 

It is for this exterior metric that we are going to construct a source. The 
motivation for this is twofold: On the one hand, it is always interesting to 
propose bounded and physically reasonable sources of gravitational fields, 
which may serve as models of compact object. On the other hand, spherical 
symmetry is a common assumption in the study of compact self-gravitating 
objects (white dwarfs, neutron stars, black holes). Therefore it is pertinent 
to ask, how do small deviations from this assumption, related to any kind of 
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perturbation (e.g. fluctuations of the stellar matter, external perturbations, 
etc.), affect the properties of the system?. However, as mentioned before, for 
sufficiently strong fields, in order to answer to this question it is necessary to 
deal with non-spherically symmetric exact solution to the Einstein equations. 

For constructing the source we shall follow a prescription given by Hernan- 
dez [11] allowing to obtain interior solutions of Weyl space-time, from known 
spherically symmetric interior solutions. Our interior solution will be ob- 
tained from the interior Schwarzschild solution (homogeneous density). 

The paper is organized as follows. In the next section we review Weyl 
solutions and the concept of relativistic multipoles. In Section 3 we describe 
the M-Q solution [10] and give its properties. The Hernandez method is 
applied in Section 4 to obtain an interior to the M-Q solution. Finally, 
results are discussed in the last Section. 

2 The Weyl metrics 

Static axysymmetric solutions to the Einstein equations are given by the 
Weyl metric [1] 



(1) 



where metric functions have to satisfy 



,zz 







(2) 
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and 



r, = A>(*?,-*y ; r,= 



(3) 



Observe that (2) is just the Laplace equation for \1> (in the Euclidean 
space), and furthermore it represents the integrability condition for (3), im- 
plying that for any "Newtonian" potential we have a specific Weyl metric, a 
well known result. 

The general solution of the Laplace equation (2) for the function \l>, pre- 
senting an asymptotically flat behaviour, results to be 



where r = (p 2 + z 2 ) 1 ^ 2 , cos 9 — z/r are Weyl spherical coordinates and 
P n (cos6) are Legendre Polynomyals. The coefficients a n are arbitrary real 
constants which have been named in the literature "Weyl moments", al- 
though they cannot be identified as relativistic multipole moments in spite 
of the formal similarity between expression (4) and the Newtonian potential. 

Another interesting way of writting the solution (4) was obtained by 
Erez-Rosen [12] and Quevedo [13], integrating equations (2), (3) in prolate 
spheroidal coordinates, which are defined as follows 



n=0 ' 



(4) 



x 



r + + r_ 
2a 



V = 



r + — r 
2a 
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x > 1 , -\<y<\, (5) 

where a is an arbitrary constant which will be identified later with the 
Schwarzschild's mass. Inverse relation between both families of coordinates 
is given by 

p 2 = aV-i)(i-y 2 ), 



z = axy. (6) 

The prolate coordinate x represents a radial coordinate, whereas the other 
coordinate, y represents the cosine function of the polar angle. 

In these prolate spheroidal coordinates takes the form 

oo 

* = S(-l) B+1 ft.Qn(x)P„(y), (7) 

n=0 

being Q n (y) Legendre functions of second kind and q n a set of arbitrary con- 
stants. The corresponding expression for the function T, has been obtained 
by Quevedo [13] 

A sub-family of Weyl solutions has been obtained by Gutsunaev and 
Manko [14] , [15] starting from the Schwarzchild solution as a "seed" solution. 
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Both sets of coefficients, a n and q n , characterize any Weyl metric [13]. 
Nevertheless, as mentioned before, these constants do not give us physical 
information about the metric since they do not represent the "real" multi- 
pole moments of the source. That is not the case for the relativistic multipole 
moments firstly defined by Geroch [17], and subsequently by Hansen [18] and 
Thorne [19], which, as it is known, characterize completely and uniquely, at 
least in the neighbourhood of infinity, every asymptotycally flat and station- 
ary vacuum solution [20], providing at the same time a physical description 
of the corresponding solution. 

An algorithm to calculate the Geroch multipole moments was developed 
by G. Fodor, C. Hoenselaers and Z. Perjes [21] (FHP). By applying such 
method, the resulting multipole moments of the solution are expressed in 
terms of the Weyl moments. Similar results are obtained from the Thome's 
definition, using harmonic coordinates. The structure of the obtained relation 
between coefficients a n and these relativistic moments allows to express the 
Weyl moments as a combination of the Geroch relativistic moments. 

2.1 The Monopole— Quadrupole solution, M — Q 

In this section we would like to describe the properties of a solution (here- 
after referred to as the M — Q solution [10]) which is particularly suitable for 
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the study of perturbations of the spherical symmetry. The main argument 
to support this statement is based on the fact that previously known Weyl 
metrics (e.g. Gutsunaev-Manko [14], Manko [15], gamma metric [8], Cur- 
zon [22], etc.) have a multipolar structure (in the Geroch sense) such that 
all the moments higher than the quadrupole, are of the same order as the 
quadrupole. In fact for the above mentioned metrics we have (odd moments 
are of course vanishing) 

Mq M = Mq R = M, 

M 2 GM = M 2 ER = —q 2 M 3 , 
15 

M 4 GM = _3Mf R = i^ 2 M 5 , 

35 

, ,,,„ 2 8r7 2 ,~o 2 4 ,,7,194 14 , 

= U." - -—MMf R = T ^^(- + jgfc), 

where q 2 is the quadrupole parameter in the Erez-Rosen metric. For the 
gamma metric results to be 



M 


= 7 M, 






M 2 


M\ 
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-7 2 ), 




M 4 


M 5 , 
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-7 2 ) 
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M 7 
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-7 2 )(1 + 
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and finally, Curzon metric is the worst case of the mentioned metrics since 
it posseses a unique parameter which represents the mass, and all higher 
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moments are proportional to increasing powers of that parameter, i.e., 



M = -a , 
M 2 = ^a 3 , 



19 - 

M 4 = an, 

105 °' 



389 7 
M 6 = al. 

6 3465 

In [10] it was shown that it is possible to find a solution of the Weyl family, 
by a convenient choice of coefficients a n , such that the resulting solution 
possesses only monopole and quadrupole moments (in the Geroch sense) 
[17]. The obtained solution (M — Q) may be written as a finite series of 
Gutsunayev-Manko [14] and Erez-Rosen [12] solutions, as follows: 

oo 

^m-q = + qV q i + q 2 m q 2 + ... = J2 <f* q <* , (8) 
where the zeroth order corresponde to the Schwarzschild solution. 

oo \2n+l 

^» = -E^— T ft»(cos0) , (9) 

n=0 ZH "T 1 

with A = Mj r and it appears that each power in q adds a quadrupole correc- 
tion to the spherically symmetric solution. Now, it should be observed that 
due to the linearity of Laplace equation, these corrections give rise to a series 
of exact solutions. In other words, the power series of q may be cut at any 
order, and the partial summatory, up to that order, gives an exact solution 
representing a quadrupolar correction to the Schwarzschild solution. 
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The simplest way to interpret physically the exact solutions obtained 
from the quadrupolar corrections described above, consists in analyzying 
the corresponding multipolar structure. Thus , it can be shown that cutting 
solution (8) at some order a, one obtains an exact solution with the following 
properties: 

• Both, the monopole and the quadrupole moments are non-vanishing: 
M = M. M-> = qM 3 . 

• All the remaining moments until order 2 (a + 1) (included) vanish. 

• All moments above the 2(a + l)-pole are of order q a+1 . Therefore, 
the solution represents a quadrupolar correction to the Schwarzschild 
solution, which is an exact solution up to the given order. 

To illustrate further our point, let us present the explicit solution up to 
first order, describing a quadrupolar correction to the monopole. The corre- 
spondig metric functions read (note a misprint in the equation (13) in [23]) 



(i) 

M-Q 




X 



3x 2 - 1 



3y 2 -l 



1 



) 




(x 2 -y 2 )(3y 2 - 1) 2 



2x 3x 



(10) 
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r (l) 

1 M—Q 



= r„o + qT„i + q 2 T Q 



+ 



225 
~24 



q 2 In 



x 2 — y 2 



225 
1024 

^(l-2/ 2 ) 



<zV - 1)(1 - y 2 )(x 2 + y 2 - 9x 2 y 2 - 1) In 2 



75 
16 
75 



16 — y 



l-^q(x 2 + 4y 2 -9x 2 y 2 + 4) 

- -a{% 2 + y 2 ) 



_^^2^ 7/ 2 W i y 2 

2 _ ,,2 4 



192 
15 



v i-y 

g 2 (2x 6 - x 4 + 3x 4 y 2 - Qx 2 y 2 + 4x 2 y 4 - y A - y 6 ) 



(x 2 — y 2 y 



1-y 2 



8 



x In 



■qx(l - y 



x — 1 



1 - W + 7y 2 - 9x 2 y 2 + 1 



x + 1 



(x 2 - y 2 ) A 
8 i 2 + 1 \ 
3 x 2 — y 2 J 

(11) 



The first twelve Geroch moments of this solution are (odd moments vanish 
because of the reflexion symmetry) 

M = M, \I> M 3 g, M 4 = 0, M 6 = -^M 7 q 2 , 



M 8 

Mia 



1060 ,, „ 2 40 ,,q o 

MV M 9 g 3 , Mio = 

3003 H 143 y ' 138567 



19880 „, n 2 146500 „, n 3 
M n g 2 H M n q 3 , 



323323 



23600 2 517600 „, 1? , 4259400,,,, 4 
-M 13 g 2 H =M 13 g 3 + M 13 g 4 . 



(12) 



437437 ' * 1062347 " * ' 7436429" 
From the expressions above, it is apparent that the parameter g = M 2 /M 3 
representing the quadrupole moment, enters into the multipole moments 
M 2n , for n > 2, only at order 2 or higher. Accordingly, solution (10)-(11) for 
an small value of q, up to order q, may be interpreted as the gravitational 
field outside a quasi-spherical source. The spacetime being represented by 
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a quadrupole correction to the monopole (Schwarzschild) solution. This is 
in contrast with other previuosly mentioned solutions of Weyl family, where 
all moments higher than the quadrupole are of the same order in q as the 
quadrupole, and therefore for small values of q they cannot be interpreted as 
a quadrupole perturbation of spherical symmetry. 

Instead of cylindrical corrdinates (p, z) , it will be useful for the next 
section to work with Erez-Rosen coordinates (r, 9) given by: 

z = (r-M) cos 9, 

p = (r 2 - 2Mr) 1/2 sinfl, (13) 
and related to prolate coordinates, by 

r 

x = T7 ~ 1> 

M 

y = cos 9. (14) 

The metric functions of the solution, up to the first order in q, hereafter 
refered as M-Q^ are: 

*ff-Q = iln(l-|)+^?(3!/ 2 -l)(3fl 2 -6fl + 2)ln(l-|) 

- (!-!)- kdfe + T^ 2 ~ w ~ ^ 
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1 r (r- i) 2 - 1 1 _ is 

2 n [(R~l) 2 -y 2 \ ~ Y 



q{l-y 2 ){R-l)\n (l 
\ (R-l) 2 + y 2 1 



1? 



2 



) 



- l) 2 - y 2 ) 2 



(16) 



with R = r/M . An study of the geodesies in this spacetime, has been recently 
presented [24] 

In the next section we shall construct a source for M-Q 1 ^ solution. 

3 An interior M-Q^ metric 

The Hernandez method [11] is based on a heuristic procedure which allows, 
starting with some spherically symmetric "seed" source, to obtain interior 
solutions describing axialsymmetric static sources, which match smoothly 
on the boundary surface, to a given metric of the Weyl family (see [11] for 
details). Some applications of this method may be found in [25]. 

For our exterior spacetime we shall choose the M — solution de- 
scribed above (up to first order in q) and our "seed" interior fluid will be the 
incompressible Schwarzschild interior solution. 

Thus, our exterior metric in Erez-Rosen coordinates read: 
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gee = e 2^ - 2 V(l-2A + A 2 sin 2 ^)=rV^ 1 - , 
9H> = e- 2 Vsin 2 ^(l-2A) = r 2 sin 2 ^ e - 2( 'V, 
gtt = -e 2 * = -(l-2M/r)e 2 ^, 
where the metric functions r„i and ty q i are given by 



V = ^(3y 2 -l)(3i? 2 -6i? + 2)ln(l-| 



(17) 



R-l 



15 



4(i?-l) 2 -y 2 16 
2 



(3y 2 -l)(i?-l) ,(18) 



"f (l-y 2 )(i?-l)ln(l 



R 

{R-lf + y 2 



(19) 



.((i?-l) 2 -y 2 ) 2 _ 

Now, Darmois conditions, in these coordinates, imply that metric compo- 
nents as well as g$e, r , gtt,r, g<j><j>,r are continuosus across the boundary surface 
(but allows a jump in g rr ^ r ). 

Thus, in the example given by Hernandez [11], the following substitutions 
on the chosen exterior metric were applied, in order to obtain the interior 
metric: 

2M/r -> r 2 /B 2 , (20) 

for the g rr metric component and the following one for the other metric 
components 



2M/r -> 1 - 



3 / rl \ 1/2 1 



.2 \ V2' 



1 - 



B 2 



2 I 1 B 2 



(21) 
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where r = is the equation of the boundary surface of the source and B 2 = 
3/(87rp s ), with p s denoting the energy density of the spherically symmetric 
"seed" solution. Since (20) does not have a continuous derivative at the 
surface, but (21) has, it is clear that Darmois conditions will be satisfied by 
the so obtained metric, on the boundary surface. 

However, in our case these substitutions lead to a metric whose g rr com- 
ponent is singular at the origin. Thus as a first step, we are going to modify 
the Hernandez rule, concerning the g rr component. In order to have a regular 
metric at the origin, we shall use (21) in those terms of g„ which produce 
the singularity at the origin and (20) in the remaining terms. For the others 
metric components we apply (21). In addition, we modify the spherical fac- 
tor in g ge and g^ by substituing r 2 with r 2 + q(r — r 2 ) 2 , in order to ensure 
regularity at the origin (however, as we shall see below this is not enough to 
assure the correct physical behaviour at the centre). Thus we obtain, 




(22) 




(23) 



9ee = (r 2 + q(r - r s ) 2 )X(r)- c ^e d ^, 



(24) 



g H = (r 2 + q(r - r s ) 2 ) sin 2 6X(r)- a{r ^e- b ^ y \ 



(25) 
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with 



X(r) = 



- ( 1 — — I 

2 I B 2 



2 \ 1/2 



2 I B 2 



2 \ 1/2' 



5 1 + X(r 



_ 15 -(l+X(r) 2 )(l-y 2 )+4X(rV 

(l-X(r)) 2 



c(r,y) 



2*1 -X(r) 
15 



4 (3y2 X) (l + X(r)) 2 -y 2 (l-X(r)) 2 
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5 



1 -X(r) 



(l-X(r))' 



+ (i-, 2 )(i-x W ) 2 a+w+y-a -*(0) a 



[(l+X(r)) 2 -y 2 (l-X(r)) 2 ] 



21 2 



1 -X(r) 



(26) 

(27) 
(28) 
(29) 



(30) 



(l+X(r)) 2 -y 2 (l-X(r)) 2 _ 

Next, let us introduce the following dimenssionles variables: a = 2M/r s and 
(3 = r/r-z, in terms of which we may write 



X{0) = 



(1-«) 1/J - (1-^) 



1/2 



n 2 



(31) 



We can now calculate the components of the energy-momentum tensor. 
However when this is done, negative values of the energy density close to 
the center appear, even in the weak field limit. To solve this problem and to 
assure a correct physical behaviour of all physical variables, we shall consider 
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the parameter q as a function of r, such that 



q = 0; (3e[0,(3 o } 



(32) 



and 



q = QoiP - Po)\P - &) 4 ; P e \J3 , (& + &)/2], 



(33) 



where (3 , (3 S and q Q are constants such that [3 D + [3 S = 2 and the value of q at 
the boundary surface, which is q [(f3 — j3 s )/2} 8 coincides with the quadrupole 
parameter of the exterior M-Q^ solution, i.e. q = q%, where q^ denotes 
the quadrupole parameter of the exterior M-Q^. It should be observed 
that since both q and its first derivative are continuous across the boundary 
surface, junction conditions are satisfied, after the replacement above. The 
specific form of (5 as well as the values of different parameters, are indicated 
in figure 1 and in its legend. 

Thus our source consists of a spherical inner core continuously matched to 
a non-spherical distribution of matter, producing a M-Q^ spacetime at the 
outside, and satisfying the continuity of the first and the second fundamental 
forms at the boundary surface. 

Calculations of different physical variables show their correct physical 
behaviour within the matter distribution. 

Expressions are extremely lengthy and therefore we omit them here, how- 
ever they are available upon request to W. Barreto. 
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The non-vanishing components are Tq, T±, T 2 2 , T|, T\. For a < 0.7 all 
these components are regular within the fluid distribution and energy condi- 
tions are satisfied (e.g. energy density is positive and larger than stresses). 
In the weak field limit (a C 1) we obtain a quadrupole correction to the 
incompressible fluid. 

4 Conclusions 

We have seen that the M-Q^ solution satisfies the requested condition to 
be considered as a quadrupole correction to the spherical symmetry, namely: 
all relativistic moments higher than the quadrupole are of higher order in q. 
Accordingly it represents, among the known members of the Weyl family, a 
good candidate to describe small departures from spherical symmetry. 

Next, we have found a source for that space-time. The interior solution 
obtained by an application of the Hernandez algorithm, matchs smoothly to 
the M-Q« metric on the boundary of the matter distribution, is regular and 
satisfies all standard physical conditions. In the weak field approximation it 
represents a quadrupole correction to the incompressible fluid sphere. 

As we have mentioned before, a spherical core was introduced in order to 
assure acceptable physical behaviour at the centre. It consists of a sphere of 
incompressible fluid with possitive energy density, larger than pressure and 
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which matches smoothly to the outer part of the source. For the models 
presented below the spherical core represents about the 0.1% of the total 
volume of the source. 

Since we are interested in slight deviations from spherical symmetry we 
shall consider models with very small values of q. Indeed the magnitude of 
the q's in the models presented are many orders of magnitude smaller than 
the values corresponding to the earth and the sun, which in our units are 
approximately 10 15 and 10 5 respectively [9]. However for a neutron star of 
one solar mass, 10 km radius and the same eccentricity as the sun, the order 
of magnitude of q is 10~ 4 [9], as in one of the models below. 

It is also worth noticing that the critical value of alfa (~ 0.7) is smaller 
than the corresponding value in the exactly spherically symmetric case (8/9) 

We have ran a large number of models for a wide range of values of 
the parameters and, both, positive and negative g's. Below we show figures 
corresponding to two models, one with small a and the other with a value of 
a close to its critical value. In both models we considered negative values of 
q since we are primarily interested in oblate objects. 

Figures 2-5 exhibit the behaviour of physical variables for an small value 
of a (weak field limit). Despite the small value of the quadrupole parameter, 
its contribution is clearly shown in figures 2-4. In figure 5 the numerical 
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error is comparable with the quadrupole contribution, and accordingly the 
latter is somehow screened by the former. 

Figures 6-8 display the physical variables for a large value of a. In this 
case, the quadrupole correction appears sharply in Tq and T\. However it is 
neglectable in T^, T| and T|. 

Finally, figure 9 display the range of parameters a and q for which the 
solution is physically acceptable. 

It is our hope that this source as well as other interiors to some Weyl 
metrics presented elsewhere ([25] and [26]) could be used as initial configu- 
rations to describe the departure from equilibrium of very compact objects 
endowed with an small but non-vanishing quadrupole structure. 
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4 . 5 




P 



Figure 1: q (multiplied by 10 4 ) as a function of (3, for (3 Q = 0.1, (3 S = 1.9 and 
g D = 10 -3 . The maximum value of q^ is 4.3046721 x 10~ 4 which we keep for 
the exterior region. 
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Figure 7: T\ = T| = T| as a function of (3 and y, for a = 6/9. All the 
surfaces overlap for q Q = and q Q = — 10~ 4 . 
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Figure 9: Dashed region indicates the range of a and q a for which models 
satisfy T ° > and T ° > T/. 
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